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' Abstract. Simulations predict that shocks from large-scale structure formation and galactic winds 

I have reduced the fraction of baryons in the warm, photoionized phase (the Lya forest) from nearly 

. 100% in the early universe to less than 50% today. Some of the remaining baryons are predicted to 

lie in the warm-hot ionized medium (WHIM) phase at T = 10^ — 10^ K, but the quantity remains 
g3 \ a highly tunable parameter of the models. Modern UV spectrographs have provided unprecedented 

^> . access to both the Lya forest and potential WHIM tracers at z ~ 0, and several independent groups 

' have constructed large catalogs of far-UV IGM absorbers along ~ 30 AGN sight lines. There is 

\ general agreement between the surveys that the warm, photoionized phase makes up 30% of 

the baryon budget at z ^ 0. Another ^10% can be accounted for in collapsed structures (stars, 
' (— I '' galaxies, etc.). However, interpretation of the ^ 100 high-ion (OVI, etc) absorbers at z < 0.5 is 

Q^, more controversial. These species are readily created in the shocks expected to exist in the IGM, but 

they can also be created by photoionization and thus not represent WHIM material. Given several 
pieces of observational evidence and theoretical expectations, I argue that most of the observed O VI 
absorbers represent shocked gas at T ^ 10^'^ K rather than photoionized gas at T < 10^'^ K, and 
they are consequently valid tracers of the WHIM phase. Under this assumption, enriched gas at 
T — 10^^*' K can account for ^10% of the baryon budget at z < 0.5, but this value may increase 
when bias and incompleteness are taken into account and help close the gap on the 50% of the 
baryons still "missing". 
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(N INTRODUCTION 

oo . 

(3 , One of the major legacies of space-based UV instruments, in particular HST and FUSE, 

is the opening of the low-redshift intergalactic medium (IGM) to spectroscopic investi- 
gation. Perversely, far more is known about the IGM at high-redshift than in the local 
universe, simply because many of the key diagnostic lines are rest-frame UV transitions 
I only accessible to optical telescopes at z > 1.6. Indeed, before HST/FOS observations 

of low-redshift AGN, there was little interest in the low-z IGM; the observed absorber 
frequency per unit redshift d^/dz of the Lya forest lines in optical spectra drops as 
(1 -\-zY''^^ at redshifts of z = 2 — 4 (Lu et al. 1991), thus predicting that Lya systems 
at z ~ would be exceedingly rare. However, FOS revealed a significant break in the 
steep power-law slope at z ~ 1.5 (Weymann et al. 1998) due to the evolution of AGN 
ionizing radiation fields and the growth of large-scale structure. The local IGM clearly 
provides astrophysical information of interest to cosmologists, and the sensitive UV 
spectrographs of the last decade have given us unprecedented opportunities. 

Large-scale simulations (e.g., Dave et al. 1999, 2001; Cen & Ostriker 1999, 2006) 
predict that the baryon fraction f^Lya in photoionized IGM (traced by Lya forest lines) 
should drop from nearly unity at z ~ 4 to a more modest fraction at local times as AGN 
radiation and large-scale structure shocks heat and ionize the IGM. Concurrently, the 



mass in the warm-hot ionized medium (WHIM, T = 10^ — 10^ K) is predicted to rise 
from nearly zero to a significant fraction of the total baryon budget at z < 0.5. The 
actual fractions in the different temperature and pressure phases depend sensitively 
on the methodology and parameters used in the simulations. Therefore, establishing 
observational values for f^Lya and f^wHiM. as well as contributions from other phases 
at z ~ 0, are of crucial importance to our understanding of everything from large-scale 
structure formation to galactic metal feedback to the distribution of AGN throughout 
cosmic evolution. 

Measuring the photoionized phase is relatively straightforward. Despite the thinning 
of the Lya forest at low z, H I absorbers are still the strongest and most common in the 
IGM. Several large surveys have established f^Lya ~ 0.3 at z < 0.5 (e.g., Penton et al. 
2000, 2004; Lehner et al. 2007; Danforth & ShuU 2008). 

Probing the WHIM phase is somewhat more complicated. Hydrogen at these tem- 
peratures is almost entirely ionized, with no diagnostic spectral lines in any waveband. 
Instead, we must rely on less direct tracers. The first is to probe the vanishingly small 
neutral hydrogen fraction (/hi < 10^^ at T ~ 10^ K) through thermally-broadened Lya 
lines in UV spectra (e.g., Richter et al. 2004, Lehner et al. 2007). This method is fraught 
with many observational ambiguities and detection challenges (Danforth et al. 2009). 

Tracing WHIM through X-ray transitions of common, highly-ionized metal ions such 
as O VII and O VIII (which peak in abundance at T > 10^ K) shows promise, but 
presents even more observational challenges. Detections of a handful of O VII and O VIII 
absorbers in the IGM have been claimed toward several AGN (Nicastro et al. 2005, Fang 
et al. 2002, 2007), but none has yet been confirmed and much doubt exists as to their 
veracity (Kaastra et al. 2006, Rasmussen et al. 2007). Current X-ray telescopes simply 
do not have the sensitivity and spectral resolution to detect the expected IGM absorbers. 

Far-UV transitions of highly-ionized metal ions (O VI XX 1032, 1038; N V AA 1238, 
1242; N VIII XXllQ, 780, etc) offer the best current hope for WHIM measurement. 
O VI is the most promising of these with ~ 100 confirmed IGM detections to date. 
O VI reaches a peak fractional abundance under coUisional ionization equilibrium (CIE, 
Sutherland & Dopita 1993) at T ~ 300,000 K, and nothing short of soft X-ray photons 
{hv > 1 14 eV) will produce it non-thermally. However, as we will see, O VI has a number 
of ambiguities as a WHIM tracer. 

The idea of tracing WHIM with O VI was pioneered in early work on individual AGN 
sight lines (Tripp et al. 2000, Savage et al. 2002, etc). Each netted a few O VI absorbers 
over a short redshift pathlength (Az < 0.6). The largest survey to date of absorption lines 
in the low-z IGM (including O VI) was carried out by Danforth & ShuU (2008, hereafter 
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DS08): a survey of 28 AGN sight lines with both FUSE and STIS/E140M data, following 
up the smaller Danforth et al. (2005, 2006) survey of 31 FUSE sight lines (z^bs < 0.15). 
Lya absorbers were identified interactively and checked for corresponding absorption 
in seven metal ion species (O VI, N V, C IV, C III, Si IV, Si III, and Fe III). The survey 
netted ~ 650 H I absorbers (many in Lyoc, Ly/3, etc. enabling curve-of-growth analysis) 
over a redshift pathlength Az = 5.2, with 83 O VI absorbers and smaller numbers of 
detections in other metal ions. The O VI detections showed d,yV/dz = 15^2 to ^ rest- 
frame equivalent width of W1032 > 30 mA, or dJVjdz = 40+g^ to V7io32 > 10 mA. 

Two other ambitious IGM O VI surveys were published concurrently with DS08 (see 
Table 1), featuring many of the same AGN sight lines. Tripp et al. (2008) searched 16 
STIS/E140M datasets blindly for O VI and surveyed H I systems for O VI absorption, 
resulting in 51 O VI detections. Despite differences in redshift pathlength definition and 
other systematic effects, they found surprisingly compatible results: dJ^/dz = 15.6]']2 4 
for W1032 > 30 mA. Thom & Chen (2008) performed a true blind O VI survey of 16 sight 
lines and found 27 O VI absorbers for a frequencjQ dJ^/dz = 10.4 ± 2.2. 

WHIM BARYONS FROM O VI ABSORBERS 

DS08 derive the statistical fraction of low-z baryons in the WHIM phase, ^^whiMj 
relative to the critical density Pcr, by using O VI absorbers as a proxy for WHIM baryons: 

V cpcr J {0/H)Q[Z/ZQ)fom JN^^y dz / 

where /ovi is the fraction of total oxygen ions in the five-times-ionized state (assumed 
to be /ovi = 0.22, the CIE peak value) and Z/Zq is the typical metallicity of the IGM, 
taken here as 10% solar. Since O VI can be produced either in v > 200 km s^^ shocks 
or by X-ray photoionization (hv > 114 eV), I introduce the parameter Qhot in Eq. 1; 
the fraction of the observed O VI that is shock-heated and thus represents legitimate 
WHIM gas. In the analysis of DS05 and DS08, 1 assumed Qhot = 1, but recent examples 
of photoionized O VI are claimed by other groups. Undoubtedly, O VI is not wholly 
generated by either thermal (2hot = 1) or photoionized (2hot = 0) means, but I present 
arguments below that most observed O VI is thermal in nature and thus the 2hot ~ 1 
approximation is appropriate. 



Shocks are Predicted to Exist 

One consistent feature of cosmological simulations is that shocks are generated by 
material falling from voids onto filaments in the cosmic web. Typical infall velocities at 
z < 0.5 are predicted to be 100-500 km s \ strong enough to produce post-shock tem- 



Thom & Chen (2008) required a detection in both OVI lines. Their d.vV/dz value is correspondingly 
lower than both DS08 and Tripp et al. (2008) and can be taken as a lower limit. 




FIGURE 1. Shocks are predicted to exist in the IGM from infall of material onto filaments and galactic 
wind feedback. Shown above are projections from a simulation by S. Skillman et al.(priv. comm.) showing 
density (left, [M0 Mpc^^]) and temperature (right, [K]). Simulation size is 128 /i^' Mpc. 



peratures of 10^ - 10^ K (Fig. lf|. If the shocked gas is sufficiently enriched, observable 
O VI absorption should be produced. Feedback from galactic winds is important to re- 
produce the observed distribution of weak O VI absorbers (DS08) through additional 
shock heating and IGM enrichment (Cen & Fang 2006; Cen & Ostriker 2006). 



IGM Cooling Times are Long 

The IGM is diffuse enough that cooling times for hot gas will be very long (Collins, 
ShuU, & Giroux 2004): 

W = ^ ~ (60 Gyr) 5,^' . (2) 

Cooling time is approximately proportional to T^''' (indexed to T(^= 10^ K) and inversely 
to metallicity, indexed to 0. 1 Z© and overdensity 850 = p/ pa = 50 (equivalent to nfj ~ 
10^^ cm^^ at z ~ 0). Shorter cooling times are possible if higher metallicity and density 
and lower post-shock temperatures are used, but it is difficult to get cooling in less than 
a Hubble time with realistic input parameters. 

Since shocks are expected from simulations and cooling times are expected to be long, 
the observed O VI gas could arise in the shock interface where gas is being coUisionally 
heated to WHIM, rather than gas radiatively cooling back to ambient temperatures. This 
explains both the velocity proximity of O VI and lower ionization species and the narrow 
range of characteristic A^ovi values. If the peak ionization species of postshock gas is 
O VII (T ~ 10^ K), we would expect to see a much broader range of column densities 
and poor correlation with H I. 



|http : //solo . Colorado . edu/^s killman /www/ Cool_Core . html] 



Poor Multiphase Correlation 



High ions studied in the far UV (e.g. O VI, N V) are usually seen within ±50 km s~ of 
the low-ionization species (H I, C III, Si III, etc). However, they exhibit poor correlation 
in column density, detection statistics, and d^/dz distribution. DS05 and DS08 show 
that the correlation between A'ovi and A'hi is weak and that the different species cover 
vastly different ranges (1.5 and 3.5 orders of magnitude, respectively). This is not true 
when comparing ions of similar species (e.g. N V/O VI or H I/Si III) which are both 
well-correlated and with similar range. 

High- and low-ionization species exhibit very different detection statistics. Comparing 
O VI, C III, and Si HI absorbers (for which the best statistics exist at low-z), we see 
that high ions and low ions appear in the same systems only ~ 25 — 30% of the time. 
More rigorously, these absorbers display a power-law behavior in absorber frequency 
per unit redshift per column density bin: d^/dz °^ N^^ . The power law slope j8 for 
low ions and high ions are noticeably distinct: jSri = 1-73 ±0.04, jSsim = 1.8 ±0.1, 
jSovi = 2.0 ±0.1 for samples of 650, 53, and 83 absorbers, respectively. This difference 
suggests a fundamental distinction between the high- and low-ion bearing gas. 



Incompatible Ionization Parameters 

Another clue pointing to the distinct origin of high and low-ionization metal species 
comes from adjacent ion pairs such as CIII/IY and Silll/IV which exhibit well- 
correlated column densities (Fig. 2). Simple CLOUDY models give a required pho- 
toionization parameter logf/ a; —1.5 ±0.5, independent of IGM metallicity, density, 
and a reasonable range of ionizing spectral energy distribution (e.g., Telfer et al. 2002). 
N V and O VI are different elements, but they correlate well in column density (Fig. 2). 
Assuming a solar abundance ratio, the same CLOUDY models require either a pho- 
toionizing radiation field an order of magnitude stronger (log?/ = —0.4 ±0.5) or high 
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FIGURE 2. Column densities of "adjacent" ion stages Si III/IV (left), C III/IV (center) and N V/O VI 
(right) from DS08 are well-correlated, but require photoionization parameters that differ by a factor of 
^ 10 between the low and high ions. Best fits (dashed) and linear relationships (dotted) are marked. 



temperatures (T = lO^ ^iO i K). 

High and low metal ions are not usually detected in the same systems, so it should be 
no surprise that they require different ionization parameters. The photoionization solu- 
tion for N V/0 VI (logt/ ~ -0.5) requires a very diffuse IGM cloud (hh ~ 10~^ cm"^) 
given what is known of the low-z metagalactic ionizing field (e.g., ShuU et al. 1999). If 
a typical O VI absorber has A'ovi = 10^^ cm^, we can solve for the line-of-sight cloud 
dimension required to account for observed, photoionized O VI: 

, ^ovi /A^ovA /O.IZqW 10-6 w 0.2 X 

An unvirialized cloud with these dimensions would display a line width of 
^abs^o ~ few X 10^ km s^^ due to Hubble broadening alone, incompatible with the 
observed O VI linewidths (bovi) = 25 ± 15 km s^^ It is possible to get the absorber di- 
mensions down to sub-Mpc scales, but only by postulating less realistic parameters (nn, 
U, Zy. higher densities (and hence stronger ionizing fields) or higher metallicities. On 
the other hand, coUisionally ionized O VI absorbers require only £abs ~ 200 — 350 kpc 
to match the observations (e.g. Stocke et al. 2006, DS08). 



Arguments for Photoionized O VI 

Tripp et al. (2008) divide their 51 O VI absorbers into "simple" and "complex" 
systems based on velocity alignment and profile shapes of O VI and H I components. 
They argue that the simple absorbers represent gas in the same phase and use the line 
widths of both species to constrain the gas temperature. In many cases, this temperature 
is less than 10^ K, suggesting that the simple O VI systems do not represent WHIM 
baryons. 

Component alignment is necessary, but not sufficient, to prove that the species are in 
the same phase. One expects pre- and post-shock material to have a different velocity 
when looking through the plane of the shock. However, sight lines parallel to the 
shock front will pass through both pre- and post-shock gas, but the velocity vector 
will be roughly in the plane of the sky. Furthermore, supernova remnants and other 
expanding shock waves show considerably stronger emission/absorption along the edges 
(limb brightening with shock velocity transverse to the line of sight) than in the center 
where a large Av is expected. The amplified column density in edge-on shocks means 
that we will see statistically more of them than face-on shocks with large velocity 
separations between pre- and post-shock gas. Given the difficulties in generating O VI 
by photoionization, this would seem to be the more logical conclusion. 

Finally, Tripp et al. (2008) classify only 37% of their O VI systems as aligned, while 
the rest show velocity separation and are deemed multiphase. Taking all the aligned 
systems as photoionized O VI, via whatever assumptions are required to make it so, this 
still leaves 2hot ~ 0.63. If a fraction of the aligned systems are also multiphase, as I have 
argued above, Qhot increases further and the Qhot ~ 1 approximation becomes valid. 



THE LOW-z BARYON CENSUS 



Thanks to concerted efforts by several groups, ~ 50% of the baryons in the local universe 
can now be accounted for as follows: 

• Stars, Galaxies, Virialized Mass: 7-10% e.g. Salucci & Persic (1999) 

• Warm Photoionized IGM (Lyct Forest): 29±4% Penton et al. (2000, 2004), 
Lehner et al. (2007), Danforth & ShuU (2008). 

• Warm-Hot Ionized Medium (O VI- WHIM): ~10%. O VI absorbers serve as a 
proxy for metal-enriched gas at T = 10^ — 10^ K (DS08), however this is probably a 
lower limit to the total WHIM fraction as discussed below. 

• "Missing": ~50% The balance of baryons may reside in one or more reservoirs 
as yet poorly explored observationally. Possibilities include gas at T > 10^ K, metal- 
poor WHIM untraceable by metal line absorption, low column density absorbers (at any 
temperature phase), and intracluster gas. 



Biases, Uncertainties, & Future Discovery Space 

Implicit in the current census of WHIM baryons are a number of uncertainties and 
biases. Most of these serve to increase ^whim and will undoubtedly account for some 
of the missing 50% of the low-z baryons once additional tracers are measured. 

• Photoionization correction (2hot < 1): As discussed above, I believe that the large 
majority of O VI absorbers observed arise in shocked gas (2hot ~ !)• However, some 
photoionized O VI absorbers may exist. Another possibility is that of a non-equilibrium 
state in which the cooling rate is faster than the recombination rate and the O VI is 
"frozen in" to a high ionization state despite being in cool gas. Both cases would result 
in non-WHIM O VI (2hot < 1) and reduce ^^whim- 

• (Z fion) uncertainty: HwfflM scales inversely with metallicity Z and ion fraction 
/ovi (Eq. 1). The assumptions of DS08 (Z = O.IZ©, /ovi = 0.22) are both likely 
upper limits for these parameters and thus t^wHiM will increase when lower-metallicity 
systems and non-peak temperature conditions are considered. 

• Lya bias: DS08 measure O VI absorption only at the redshifts of Lya "signposts". 
Thus, the population of high-ion absorbers without Lya is unaccounted for in the present 
census. Both Tripp et al. and Thom & Chen perform "blind" surveys. Tripp et al. 
find several "naked" O VI systems, but all are plausibly intrinsic to the AGN rather 
than intervening IGM. Thom & Chen find one out of 27 O VI absorbers with no Lya 
counterpart; a close look at these data shows Lya at low significance. I conclude that the 
Lya bias is likely small and has only a few percent effect on flwHiM- 

• Metallicity bias: Low-metallicity and metal-free WHIM absorbers are likely if 
large-scale structure formation shocks are a major source of IGM heating. Void metal- 
licity is seen (Stocke et al. 2007) to be at least an order of magnitude lower than the 
canonical 10% measured in DS08. Galactic feedback shocks may account for a large 
fraction of what we see in O VI and other metals. Tracing low metallicity WHIM will 
require broad Lya lines, with all the observational difficulties they represent. 



• Temperature bias: Each FUV WHIM tracer is sensitive only to a narrow range of 
temperatures, with O VI only common over T — 10^ — 10^ K. At high T , the plummeting 
hydrogen neutral fraction and rapidly broadening line profile make broad Lya lines 
sensitive to only the coolest WHIM. Given the shape of the cooling curve, hotter WHIM 
is expected to be more common, but it remains unmeasured in the current version of the 
baryon census. 

• Weak absorbers: Finally, the observed A'ovi distribution d^/dz A'^yi "i^ans that 
equal mass is contributed from many weak absorbers as from few strong ones. Limits in 
sensitivity mean that there may be large reservoirs of WHIM in weak absorbers: A^ovi < 
10^3 cm"2 (Wmi < 12 mA). High S/N COS observations of weak O VI absorbers at 
z > 0.1 will help to quantify this uncertainty. 
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